Tricritical point of J^-J^ Ising model on hyperbolic lattice 
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A ferromagnetic-paramagnetic piiase transition of the two-dimensional frustrated Ising model 
on a hyperbolic lattice is investigated by use of the corner transfer matrix renormalization group 
method. The model contains ferromagnetic nearest-neighbor interaction Ji and the competing 
antiferromagnetic interaction J2. A mean- field like second-order phase transition is observed when 
the ratio k — J2/J1 is less than 0.203. In the region 0.203 < k < 1/4, the spontaneous magnetization 
is discontinuous at the transition temperature. Such tricritical behavior suggests that the phase 
transitions on hyperbolic lattices need not always be mean-field like. 
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I. INTRODUCTION 

The Ising model on the Cayley tree is known by its 
singular property, where the magnetic susceptibihty of 
the spin at the root of the tree diverges at a temper- 
ature T^ despite the fact there is no singularity in the 
partition function of the whole system [1] . This is a kind 
of phase transition which can be explained by the Ising 
model on the Bethe lattice. It has been known that the 
Ising model on hyperbolic lattices, which arc negatively 
curved in the two-dimensional (2D) space [3|, exhibits 
similar aspects in common [3|, y, |5|. The universality 
class of the ferromagnetic-paramagnetic phase transition 
of this model has been so far considered to be mean-field 
like. Recent numerical studies have supported this con- 
jecture |i,0,i,ii- 

In this paper we study effects of the antiferromag- 
netic next-nearest- neighbor (NNN) interaction J2 , which 
competes with the ferromagnetic nearest-neighbor (NN) 
one J^ , on the ferromagnetic-paramagnetic phase tran- 
sition of the 2D Ising model on a hyperbolic lattice. 
We use the corner transfer matrix renormalization group 
(CTMRG) method [H, [H, [3, which is a variant 
of the density matrix renormalization group (DMRG) 
method [1^, [13, H, [l3|i for the calculations of ther- 
modynamic functions. As we show in the following, 
the transition temperature Tq monotonously decreases 
with the frustration parameter k — J2/J1 in the region 
< K < 1/4, where the ground state spin configuration is 
completely ferromagnetic. We find that there is a tricrit- 
ical point when the parameter k is equal to k^ = 0.203. 
The ferromagnetic-paramagnetic phase transition is of 
the second order for < k < k^,, whereas it turns into 
the first order one for k^ < k < 1/4. 

In the next section, we explain the so-called (5, 4) lat- 
tice in the 2D hyperbolic space and introduce the Ising 
Hamiltonian on it. As a theoretical ideal, we consider 
phase transition on the Bcthc lattice with the coordina- 
tion number four, which is equivalent to the (c»,4) hy- 
perbolic lattice. In Sec. Ill we present numerical results. 




FIG. 1: The (5,4) hyperbolic lattice drawn in the Poincare 
disc. The open circles represent the Ising spins sites. The 
next-nearest-neighbor interactions are here represented by the 
dashed lines inside the pentagons. 



Temperature dependence of free energy, spontaneous and 
induced magnetizations are shown. We analyze these 
thermodynamic functions around the transition temper- 
ature Tq for several values of k, and determine the critical 
exponents a, /?, and 5. We summarize the observed phase 
transition in the last section. 



II. FRUSTRATED ISING MODEL ON 
HYPERBOLIC LATTICE 

We consider a hyperbolic 2D lattice shown in Fig. [TJ 
where four pentagons share their apexes. Such lattice 
is conventionally called as the (5,4) lattice, where the 
number five represents number of the sides of each pen- 
tagon and the number four is the coordination number. 
Consider the Ising model on this lattice, where on each 
lattice site labeled by i there is an Ising spin variable 
CTj = ±1. We assume ferromagnetic interactions between 
NN spin pairs shown by the full lines in Fig. [1] and the 
antiferromagnetic interactions between NNN pairs shown 




FIG. 2: Ground state spin configurations for k < 1/4 (left) 
and K > 1/4 (right) on the (oo,4) lattice, which coincides 
with the Bethe lattice with the coordination number four. 
Note that only a finite number of spins is here depicted from 
the (oo, 4) lattice. Open circles represent spin variable a = +1 
whereas the full circles correspond to a = — 1. 



by the dashed hnes. 
represented as 



The Hamiltonian of the system is 



H 



(jj)=NN 



{ifc)=NNN 



(1) 



where J^^ > is the ferromagnetic couphng constant be- 
tween the NN pairs {ij) and J2 > is the antiferromag- 
netic one between the NNN pairs {ik). Let us define 
a parameter k = J2I J\ that represents strength of the 
frustration. 

For purpose of obtaining brief insight of the phase 
structure of the (5, 4) Ising model, we observe the model 
from a wider framework. Let us introduce the Ising 
model on the (n, 4) lattice where four 71-gons (polygons 
of the n"^ order) meet at each lattice point. When the 
number of sides n (> 5) is multiple of four including 
the case n = 00, the ground-state spin configuration at 
zero temperature is easily obtained. Figure 2 shows the 
ground-state configurations for the case n = 00, where 
the lattice is nothing but the Bethe lattice with the co- 
ordination number four. For the complete ferromagnetic 
configuration shown on the left, the energy expectation 
value per site is 



-■ Fern 



-2 Ji -f 4J2 



(2) 



and for the 'up-up-down-down' structure shown in the 
right, the value is 



- uudd 



= -4J, 



(3) 



Therefore, the energy cross over epono = ^uudd ^^ located 
at J^ = 4J2, equivalently at k = 1/4. This ground state 
alternation is common for all the cases where n (> 5) 
is multiple of four. If not, the ground state spin config- 
uration for large k is not unique and is probably disor- 
dered. In the case of the (5, 4)-lattice, one of the ground 
states in the large k region can be constructed by joining 
the pentagons with either 'up-up-up-down-down' or 'up- 
up-down- down- down' spin configurations. After a short 



algebra, one obtains the energy per site 



_ 2 12 

'uuudd ^uuddd ~ ~^ 1 ^ 



(4) 



for the assumed configurations. Hence the energy cross 



over £p^„„ = £ 



uuudd 



also occurs at J-^ = 4 J2 



At finite temperature, the ferromagnetic-paramagnetic 
phase transition is observed in the small k region 20]. 
Consider a single-site mean-field approximation on arbi- 
trary (n, 4) lattice. A mean field variable h is expressed 
as 



h = (-4Ji + 8J2) (ct) = -(4 - 8k) Ji (a) 



(5) 



where {a) is the expectation value of the Ising spin. 
A self-consistent condition for (a) leads to the ferro- 
magnetic-paramagnetic phase transition with the criti- 



cal temperature T^ 



M.F./ 



(4 — %K)J^lk^, where A:^ 



is the Boltzmann constant. Within this approximation, 
the transition is always of the second order in the region 
< K < 1/4, since the effect of J^ appears as the rescal- 
ing of the mean-field h as in Eq. (5). It should be noted 
that T^-^\k = 1/4) = 2Ji/fcB is larger than zero. The 
mean-field approximation predicts another ordered state 
in the region k > 1/4, where the 'up-up-down-down' spin 
configuration is favored if the lattice geometry allows the 
ordering. 

An improvement to the mean-field approximation is 
achieved by increasing the number of sites that are not 
averaged. The simplest case is the Bethe approximation, 
which treats additional spins a-^^ a2, cTg , and a^ that sur- 
round the central site cr, as shown in Fig. 2 (left). On 
arbitrary (n, 4) lattice, the mean field for the surround- 
ing four spins Ui, a 2, Cg, and a^ is given by 



/i,-(-3Ji-f6J2)(a) 



(6) 



As an effect of the next-nearest-neighbor interaction, the 
central spin a also feels the mean field, however, of dif- 
ferent strength, 



h = 8^2 {(^) 



(7) 



in addition to the direct ferromagnetic interaction with 
the surrounding spins —J^a [a^ + (T2 + o's + o'4)- Con- 
sidering these interaction terms, one obtains the self- 
consistent relation 

(c^) = ^ E cr exp[-/3/i„(cri + (T2 + cTg + 0-4) - i3h^, a 
+(3Ji {a I + (72 + (7-^ + (7^) a 

-/3J2 (cri(T2 +'^2Cr3 + (73(74 + fT40-i)] , (8) 

where /3 = 1/k^T and where Z is the partition function 
Z = J2 exp[-/3/i„(cri +(72 +^^3 +(74) -/3/ibCr 

+(3Jl ((7^ + (72 + (73 + (74) (7 

-/3J2 (cti<72 + '^20-3 + (73(74 + a^<7^)] . (9) 




nr^ 



FIG. 3: Spontaneous magnetization (cr) obtained by the 
Bethe approximation. 




1.05 



FIG. 4: Spontaneous magnetization of the J1-J2 Ising model 
on the (00, 4) lattice. 



The configuration sums in Eqs. ([8]) and ([9]) are taken over 
the five spins a, cti, CT2, (73, and 174. The factorization 



VF(crj,cr) = exp 



-PKa, + (3J,a,a - (3^a 



(10) 



for i = 1, 2, 3, and 4 further simplifies the expression so 
that the partition function has the form 

Z = ^M/((Ti,a)M^(a2,a)M^(a3,a)W^(a4,a) 

exp[-/3J2 (criCr2 +<^2<^3 + '^3<^4: + '^4'^l)] -(H) 

Since it is not a trivial task to find out an analytic so- 
lution of the self-consistent Eqs. (6)-(9), we solved them 
numerically. We use the parametrization Ji = 1 and 
fcg = 1 throughout this article in the numerical calcula- 
tions. Figure 3 shows the calculated spontaneous magne- 
tization M = {a). The second-order phase transition is 
detected in the whole region < k < 1/4. As observed in 
the single-site mean-field approximation, the transition is 
of the second order, and the transition temperature re- 
mains finite even at k = 1/4. Further improvement of 
the Bethe approximation can be achieved by means of 
gradual increase of unaveraged spin sites. A series of 
such approximations is known as the coherent anomaly 
method (CAM) [2l|. Here, we do not proceed with the 
CAM analysis; we perform extensive numerical calcula- 
tions by the CTMRG method instead. 

It has been known that the spin expectation value (cr) 
can be calculated exactly at the root of the Cayley tree, 
which can be treated as the Bethe lattices [ij. For the 
frustrated Ising model on the (00, 4) lattice shown in 
Fig. 2, the expectation value is expressed as 

<^) = ^^aW^'(ai,a)I¥'(a2,a)I^'(a3,a)I^'(a4,a) 

exp[-/3J2(criCr2 + 0-20-3 + '^S'^A + '^4<^l)] (12) 

with the definition of the effective partition function 



Z' 






W'{a^,a)W'{a^,a)W'ia3,a)W'{a^, a) 



exp[-/5J2(criCT2 + (72^3 + cr3Cr4 + a^a^^)] , (13) 

where the new factor W'{a^,a) represents a Boltzmann 
weight for a branch which connects the root spin a with 
the nearest spin site cji. (c.f. Fig. 2.) This new factor 
W'{a^,a) can be calculated from W{a^,a) in Eq. (10) 
repeating the application of the recursive transformation 

exp[/3Ji cr.cr - /3J2 (crsi + s^s^ + s^s^ + s^a)] 

(14) 

for many times until it converges [l|. Thus W'{a^,a) 
contains the effect of distant sites on the Bethe lattice. 
Figure 4 shows the spontaneous magnetization AI — (a) 
calculated by Eq. (12) using W'{<7^,(t) numerically ob- 
tained from Eq. (14). The transition is of the second 
order in the region < k < 0.183 and is of the first or- 
der in 0.184 < K < 1/4. One can carry out perturbative 
calculations to ensure that the transition temperature 
on the Bethe lattice is zero a.t k — 1/4. The difference 
between Fig. (3) and Fig. (4) comes from the effect of 
distant interacting spin sites, which might be essential in 
the tricritical behavior on the (5, 4) lattice as studied in 
the next section. 



III. NUMERICAL RESULTS BY CTMRG 

In this section we analyze the thermodynamic prop- 
erty of the Ising model on the (5, 4) lattice by use of the 
CTMRG method [H, [M |T1. The method is a variant 
of the DMRG method [16|, IH ^ appHed to 2D classical 
models [22l |. It has been known that the partition func- 
tion Z of a square-shaped finite-size system can be calcu- 
lated as a trace of the fourth power of the so-called corner 
transfer matrix (CTM), which represents a Boltzmann 



Ai 
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0.05 



0.25 



K = Jj^Ji 



FIG. 5: Phase diagram of the J1-J2 Ising model on the (5,4) 
hyperbolic lattice. 



weight of a quadrant of the whole system [1| . Although 
the matrix dimension of the CTM, which is denoted by 
C, increases exponentially with the linear size of the sys- 
tem, it is possible to transform it into a renormalized one 
C with a smaller matrix dimension m [23[ by means of 
the RG transformation obtained from the diagonalization 
of p = C"' or C [l|, [ij, [1^ . This transformation is not ex- 
act but is highly accurate in such sense that Z = TrC* 
is a good approximation oi Z — Tr C^ . One can pre- 
cisely calculate thermodynamic (or one-point) functions, 
such as the free energy F — —k^T log Z and the spon- 
taneous magnetization M, for a sufficiently large finite- 
size system by use of the CTMRG method. Since the 
(n, 4) lattice can be divided into four equivalent parts 
(the quadrants), which share the central site a on their 
edges, it is also possible to apply the CTMRG method 
to statistical models on these lattices [g, |^, [I4I . 

In order to study critical phenomena correctly on hy- 
perbolic lattices, we put the following remarks. We al- 
ways consider a latices system whose linear size L is 
several times larger than the corresponding correlation 
length ^, so that the central site a is sufficiently away 
from the system boundary. The lattice sites in the 
area within the distance of the order of ^ from the sys- 
tem boundary are affected by the imposed ferromagnetic 
boundary condition, where all the Ising spins at the sys- 
tem boundary point to the same direction. It should 
be also noted that portion of such sites that are 'near 
the boundary' in the hyperbolic geometry remains finite 
even in the thermodynamic limit L — > c» jlOl . Illl | , where 
the situation is similar to the case of the Cayley tree [l| . 
Disregarding all these sites 'near the boundary', we focus 
on the thermodynamic properties of the Ising spins deep 
inside the system [3|, [1^ . 

Figure [5] shows the phase diagram of the system in 
the parameter region < k < 1/4. The ferromagnetic- 
paramagnetic phase boundary is determined from the 
temperature dependence of the free energy F(k; T) and 
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FIG. 6: Dependence of the free energy F(k; T) on tempera- 
ture when (a) k = 0.18 and (b) k = 0.22. 




FIG. 7: Spontaneous magnetization M for k = 0. 

the spontaneous magnetization M{k;T) which we show 
in the following. As an effect of the competing inter- 
actions, the transition temperature Tq (k) monotonously 
decreases with k towards rQ(l/4) = 0. In the region 
< K < 0.203 the transition is of the second order. In 
contrast, when 0.203 < k < 1/4, we observe a first-order 
transition; the tricritical point is located at k^ = 0.203. 
Figure [HI shows the free energy F{k;T) at k = 0.18 and 
K — 0.22. In the region 0.203 < k < 1/4, the free energy 
F{k:T) is not a differentiable function at the transition 
temperature, as shown in Fig. 6(b). 

Figure [7] shows the spontaneous magnetization M and 
its square, M^, when there is no frustration k — 0. The 
squared magnetization M^ is proportional to Tq — T, the 
behavior which agrees with the critical exponent f3 = 1/2. 
In the second-order transition region < k < 0.203, the 
universality class remains mean-field like and the magne- 
tization curve satisfies the scaling form 

MiK;T)^B{K)[mK)~Tf^ (15) 

around the transition temperature Tq{k). The prefactor 
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FIG. 8: Inverse of the prefactor B{k), which characterizes the 
mean-field like transition observed in the spontaneous mag- 
netization. 



FIG. 10: The discontinuity function of the spontaneous mag- 
netization at the transition temperature Tq in the first-order 
transition region k^ < k < 1/4. 
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FIG. 9: Temperature dependence of the spontaneous magne- 
tization M{hi]T) for several values k around n^. 



B{k) is an increasing function of k, and diverges at certain 
point K = K(.. Figure [S] shows the inverse of the prefactor 
B{n) which linearly decreases to zero in the vicinity of 
K^. We obtain k,^ — 0.2027 from the linear fitting. 

One can also estimate k^ out of the discontinuity in 
the spontaneous magnetization M{k;T) in the region 
0.203 < K < 1/4. Figure [9] shows M(k; T) around n = n^. 
We calculate the discontinuity function (or the jump in 
the magnetization) D{k) = M{k;T^)^ where Tq~ cor- 
responds to temperature just below the transition tem- 
perature Tq{k). As shown in Fig. [TOl the discontinuity 
function satisfies relation 



D{k) oc (k - kJI/^ 



(16) 



around k = k^,. Performing a linear fitting shown by 
the dashed lines, we obtain k^ — 0.2033. Comparing 
this value with k^ = 0.2027 obtained from the data in 
second order region, we conclude that the tricritical point 
is located at k, = 0.203 ± 0.001. 



The observed tricritical behavior around k — k^ in in 
accordance with the Landau free energy 



F(M, t)^a M^ + b{K^- k)M^ + ct M^ 



(17) 



where a, b, and c are positive constants or slowly varying 
functions of temperature. In the second order transition 
region k < k^, the second and the third terms in F{M, t) 
are dominant in the vicinity of the phase transition, and 
the parameter t coincides with [T — Tq{k)]/Tq{k). Ne- 
glecting the first term in F{M,t) below T < Tq{k), we 
can obtain the spontaneous magnetization M that mini- 
mizes F{M, t) from the equation 46 {k^ — k)M'^ -\- 2ct = 0. 



The behavior M'^ 



It I coincides with the numerical re- 



sult shown in Fig. 7. In the first order transition region 
K > Kj., all three terms in F(M, t) are important for the 
minimum of the free energy. After short calculations, 
one can confirm that the jump of the spontaneous mag- 
netization at the transition temperature coincides with 
Eq. (16), which we have verified from the numerical data 
shown in Fig. 10. 

At the tricritical point k — k^, the second term in 
F{M,t) in Eq. (17) vanishes. Thus, the spontaneous 
magnetization is determined from 6aM'^ + 2ct = 0, and 
M^ is proportional to T — Tq. The dependence of M^ 
with respect to temperature T obtained from the numer- 
ical calculation is shown in Fig. [TTJ The data are in 
accordance with Af oc (T — Tq)^/^, which corresponds 
to the exponent /3 = 1/4 at tricriticality. In the same 
manner it is expected that the specific heat diverges as 
(Tq — r)^^/^, which corresponds to the critical exponent 
a = 1/2. Figure [T2] shows the numerically calculated spe- 
cific heat at k^ — 0.203, which agrees with the expected 
temperature dependence. For comparison, in the inset of 
Fig. [m we show the data at k = 0.18, which agrees with 
a = 0. 

Effect of an external magnetic field H may be included 
in the Landau free energy by adding the interaction term 




FIG. 11: The fourth power of the spontaneous magnetization 
around k = k^. 
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FIG. 13: Induced magnetization (i) at the tricritical point 
n — He and (ii) inset: that at the transition temperature 
when K < K„. 
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FIG. 12: Specific heat around the tricritical point. 



dHM to F(M,t) in Eq. (17), where d is a positive con- 
stant. From the assumed form of the free energy it is ex- 
pected that M^ and M^ are, respectively, proportional 



to Tn — r when k < k^ and when n — 



For the 



confirmation, we observe the induced magnetization at 
criticahty when k < k^. Figure [T3] shows the induced 
M with respect to H. In the region of the second or- 
der phase transition, we obtained the magnetic exponent 
(5 = 3 as expected. However, the value of the exponent 6 
is around 7 at the tricritical point, not 5 as expected from 
the Landau free energy F{M, t) + dHM. Such pathologi- 
cal behavior in the induced magnetization is a remaining 
piece of the puzzle of the current study on the (5, 4)- 
hyperbolic lattice. Speak about the J^ ~ J2 Ising model 
on the (oo,4)-Bethe lattice, we obtain S ^ 6 a.t the tri- 
critical point from numerical calculations. Future studies 
on {n, 4)-lattices for n > 6 would provide a hint to these 
unexpected values of d. 



IV. CONCLUSIONS 

We have studied a ferromagnetic-paramagnetic phase 
transition of the J1-J2 Ising model on the (5, 4) hyper- 
bolic lattice. A tricritical point has been found when the 
ratio K = J2/ Ji, which represents the strength of frustra- 
tion, is equal to 0.203. It should be noted that the pres- 
ence of the first-order transition cannot be obtained by 
the single-site mean-field approximation applied to this 
system. This is in contrast to the known fact that the 
phase transition of the nearest-neighbor IsingjPotts, and 
clock models exhibits mean- field nature [3, [a, 0, la, M, \lM ■ 

The observed second order phase transition in the pa- 
rameter region k < k^ = 0.203 belongs to the mean-field 
universality class, which is characterized by the expo- 
nents a = 0, /3 = 1/2, and 5 = 3. At the tricritical point 
we observe a = 1/2, P = 1/4, which are in accordance 
with the Landau free energy written in the even polyno- 
mial of order parameter. The observed value of exponent 
(5 « 7 at the tricritical point is the only exception which 
requires further detailed studies. 

As an effect of the frustration, entropy of the ordered 
phase shall be enhanced compared with the ordered state 
that has the same spontaneous magnetization under J2 = 
0. We conjecture that this enhancement effect creates a 
minimum in the Landau free energy, which may be the 
reason of the first order transition we have observed here. 
The tricritical point is also present in the (00, 4) lattice, 
which is nothing but the Bethe lattice. This suggests that 
the suppression of the loop back effect in the hyperbolic 
lattice is essential for the appearance of the tricritical 
behavior. 

Determination of the phase diagram in the region 
K > 1/4 is challenging because the ground state spin 
configuration becomes non-trivial as has been discussed 
in Section II. Because the (5,4) hyperbolic lattice con- 
sists of pentagons, the lattice does not decouple into sub- 
lattices even when J^ — 0. For the study of this region, 
we have to modify the CTMRG algorithm in order to 



treat ordered states with non-trivial spin patterns. 
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